Multi-level Monte Carlo finite volume method for shallow water equations with uncertain parameters applied to landslides-generated tsunamis
نویسندگان
چکیده
Two layer Savage-Hutter type shallow water PDEs model flows such as tsunamis generated by rockslides. On account of heterogeneities in the composition of the granular matter, these models contain uncertain parameters like the ratio of densities of layers, Coulomb and interlayer friction. These parameters are modeled statistically and quantifying the resulting solution uncertainty (UQ) is a crucial task in geophysics. We propose a novel paradigm for UQ that combines the recently developed IFCP spatial discretizations with the recently developed Multi-level Monte Carlo (MLMC) statistical sampling method and provides a fast, accurate and computationally efficient framework to compute statistical quantities of interest. Numerical experiments, including realistic simulations of the Lituya Bay mega tsunami, are presented to illustrate the robustness of the proposed UQ algorithm.
منابع مشابه
Uncertainty quantification in tsunami modeling using multi-level Monte Carlo finite volume method
Shallow-water type models are commonly used in tsunami simulations. These models contain uncertain parameters like the ratio of densities of layers, friction coefficient, fault deformation, etc. These parameters are modeled statistically and quantifying the resulting solution uncertainty (UQ) is a crucial task in geophysics. We propose a paradigm for UQ that combines the recently developed path...
متن کاملMultilevel Monte Carlo Finite Volume Methods for Shallow Water Equations with Uncertain Topography in Multi-dimensions
The initial data and bottom topography, used as inputs in shallow water models, are prone to uncertainty due to measurement errors. We model this uncertainty statistically in terms of random shallow water equations. We extend the Multi-Level Monte Carlo (MLMC) algorithm to numerically approximate the random shallow water equations efficiently. The MLMC algorithm is suitably modified to deal wit...
متن کاملMulti-level Monte Carlo finite volume methods for nonlinear systems of conservation laws in multi-dimensions
We extend the Multi-Level Monte Carlo (MLMC) algorithm of [19] in order to quantify uncertainty in the solutions of multi-dimensional hyperbolic systems of conservation laws with uncertain initial data. The algorithm is presented and several issues arising in the massively parallel numerical implementation are addressed. In particular, we present a novel load balancing procedure that ensures sc...
متن کاملA Coupled Numerical Model for Tsunamis Generated by Subaerial and Submarine Mass Failures
This paper presents a new numerical model simulating tsunamis generation by landslides. Water, air, and the slide (considered either as a viscous fluid or as a rigid bloc using the penalization method) are described by Navier-Stokes equations, expressed in a unified single fluid approach. The PLIC-VOF method is used to describe the motion of fluid interfaces. We present results for two test cas...
متن کاملSimulation of the Propagation of Tsunamis in Coastal Regions by a Two-Dimensional Non-Hydrostatic Shallow Water Solver
During the last years, great effort has been addressed by several authors to simulate the propagation of solitary waves/ tsunamis, tides or surges, due to the tremendous damages and losses of human lives in the inundated rural and residential areas. Tsunamis are sea waves usually generated by undersea landslides and earthquakes. They can be regarded as long/solitary waves with small amplitude a...
متن کامل